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FIGURE  CAPTIONS 


(1)  The  modulus  of  the  nine  surface  modes  (slopes)  i? 
plotted  at  five  different  times,  i.e.,  five  "snap¬ 
shots"  of  the  surface  spectrum. 

(2)  The  modulus  of  the  nine  surface  modes  (slopes)  are 
plotted  as  continuous  functions  of  time.  At  any 
instant  of  time  the  surface  spectrum  can  be  read  from 
the  labeled  curves. 

(3)  The  surface  distortion  for  the  tank  experiment  is 

plotted  and  compared  with  three  calculations.  Two 
calculations  are  at  the  resonance  c  =  c^  ,  (i)  linear, 

(ii)  modal,  and  the  third  calculation  is  modal  theory 

off  resonance. 

(4)  The  continuous  growth  of  the  surface  spectrum  is  indi¬ 
cated  for  the  calculation  depicted  in  Figure  (2),  but 
with  the  four  surface  wave  interaction  "turned  on". 

(5)  The  continuous  growth  of  the  surface  spectrum  is  indi¬ 
cated  for  the  calculation  depicted  in  Figure  (4) ,  but 

with  additional  modes  situated  at  k  t  K/2  and  k  ±  4K  • 

o  o  2 

(6)  The  development  of  the  ocean  surface  spectrum  is  indi¬ 
cated  at  four  distinct  times.  The  discrete  nature  of 
the  calculated  spectrum  is  approximated  by  the  continuous 
line  joining  the  21  modes. 

(7)  The  same  as  Figure  (6)  but  at  latter  times. 

(8)  The  surface  listortion  of  the  ocean  surface  is  indi¬ 
cated  by  the  straight  line,  i.e.,  no  distortion. 
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1.  INTRODUCTION 


The  first  report  in  this  series  [Watson,  West  and  Thomson,  (I), 
PD-72-030  (RADC-TR-73-74) ]  was  concerned  with  the  development 

of  a  formalism  which  would  allow  a  direct  calculation  of  the 
non-linear  interaction  among  surface  gravity  waves.  The  tech¬ 
nique  used  was  to  expand  the  ocean  surface  in  a  series  of 
modes,  i.e.,  a  finite  Fourier  series,  and  express  Bernoulli’s 
equation  and  the  kinematic  boundary  condition  at  the  ocean 
surface  in  terms  of  interacting  modes.  These  modal  equations 
formed  a  system  of  rate  equations  with  quartic  and  triplet 
mode  couplings  giving  the  non-linear  interactions.  These 
equations  were  then  transformed  to  normal  coordinates,  i.e.,  to 

the  eigenmodes  of  the  linear  system,  and  integrated  numerically 
on  a  computer.  A  calculation1'2  with  simulated  initial  conditions 
compared  favorably  with  the  Ben j amin-Feir  (BF)  experiment  (1967). 

In  this  second  report  we  wish  to  extend  the  above  analy¬ 
sis  to  include  the  interaction  of  the  surface  gravity  waves 
with  a  surface  current.  We  will  treat  the  surface  current  as 
being  generated  by  a  "massive"  internal  wave,  i.e.,  an  internal 
wave  which  does  not  react  to  the  surface  gravity  waves  [see 
Thomson  and  West,  PD-72-023  (RADC-TR-72-280) ] ,  and  consider 
only  the  linear  interaction  between  the  current  and  surface 
gravity  waves.  The  consequences  of  relaxing  this  assumption, 


i.e.,  allowing  the  internal  wave  to  change,  is  being  inves¬ 
tigated  elsewhere  (Watson,  West,  Cohen,  PD-73-032)  .  Section  2 
reviews  the  construction  of  the  linear  equation  and  obtains 
results  consistent  with  those  of  M.  Milder  (1972)  and 
indicates  an  error  by  Zachariasen  (1972)  pertaining  to  the 
surface  on  which  the  equations  are  evaluated. 

In  Section  3,  a  translating  surface  current  is  expanded 
in  harmonics  and  introduced  into  the  linear  equations,  which 
are  expressed  in  terms  of  modes.  The  equations  are  written 
in  two  dimensions,  since  the  added  complexity  over  the  one 
dimensional  equations  is  not  great.  The  eigenmode  equations 
close  to  resonance  are  rather  simple,  a  given  mode  £  being 
coupled  only  to  its  neighboring  modes  where  £  is  the 
fundamental  wavenumber  of  the  internal  wave.  A  hermitian 
representation  of  the  interaction  equation  can  be  approxi¬ 
mated  near  resonance  and  an  interesting  form  of  tnis  equation 
is  derived  in  two  dimensions. 

The  two  dimensional  nature  of  the  equations  is  not 
exploited  in  this  report,  but  a  discussion  indicating  the 
possible  direction  of  future  work  is  given.  In  the  one 
dimensional  form,  the  expressions  are  compared  to  those 
developed  by  Rosenbluth  (1971)  ,  where  an  inconsistency  in 
the  level  of  approximation  is  found.  This  inconsistency 
arises  from  the  additional  terms  uncovered  in  the  analysis 
of  Zacharisen  (1972).  A  modal  expression  exact  to  order  1/N 
is  derived  and  the  solution  in  terms  of  convolutions  with 
discrete  Greens  functions  is  formulated. 
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The  modal  equations  are  solved  numerically  in 
Section  4  for  a  number  of  cases.  The  first  calculation 
employs  initial  conditions  which  simulate  the  experiment  of 
Lewis  (1971).  In  this  experiment,  an  internal  wave  was 
generated  at  a  fluid  interface  near  the  bottom  of  a  tank. 

An  infinitesimal  surface  wave  was  mechanically  generated 
and  the  modulation  of  the  surface  height  and  slope  recorded. 
We  calculate  the  modulation  as  a  function  of  time,  as  well 
as  the  development  of  the  surface  spectrum  and  compare  the 
results  with  the  experimental  data, 

A  second  calculation  is  presented  in  which  the  para¬ 
meters  are  characteristic  of  the  ocean  environment  rather 
than  those  of  a  tank.  The  case  of  a  single  (dominant)  sur¬ 
face  wave  in  an  unsaturated  surface  spectrum  interacting 
with  an  internal  wave  is  considered.  The  results  are  dis¬ 
cussed  in  the  context  of  the  linear  theory  developed  in 
Section  3  and  the  effect  of  nonlinear  tranfer  of  energy 
between  surface  waves  is  also  considered.  The  distinction 
between  the  interaction  regimes  for  tank  experiments  and 
the  ocean  experiments  is  also  discussed.  The  important 
quantity  is  determined  to  the  relative  magnitude  of  the 
correction  to  the  dispersion  relation  for  the  eigenmodes 
and  the  coupling  coefficients  of  the  surface  waves  to  in¬ 
ternal  waves. 
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2.  LINEAR  WAVE  EQUATION 

This  section  is  essentially  a  review  of  the  derivation 
of  the  linear  wave  equation  and  may  be  omitted  by  those  fa¬ 
miliar  w:th  its  construction.  In  this  development  we  consider 
the  surface  gravity  waves  to  be  infinitesimal  and  to  be  inter¬ 
acting  with  a  non-uniform  surface  current.  The  surface  cur¬ 
rent  may  be  generated  by  a  non-dispersive  internal  wave,  a 
swell  passing  through  the  region  of  interest  or  by  a  number 
of  other  mechanisms.  For  the  present  we  will  not  be  concerned 
with  how  the  surface  current  is  produced,  only  that  it  is 
present.  This  analysis  compliments  that  in  PD  72-023  (RADC- 
TR-72-280) ,  i.e.,  the  linear  wave  equation  being  constructed 
includes  terms  which  were  considered  negligible  previously  and 
found  by  Zachariasen  (1972)  to  be  important. 

We  assume  the  ocean  to  be  both  homogeneous  and  irrota- 
tional  so  that  we  can  use  a  potential  description  of  the  velo¬ 
city  field,  i.e.,  since  Vxu  =  0  we  can  write  u  =  V4>.  The  ocean 
is  also  assumed  to  be  incompressible  (V*u  =  0)  so  that  the 
velocity  potential  satisfies  Laplace's  equation, 

V2*  =  0  .  (2.1) 

The  integral  of  the  momentum  equation  yields  Bernoulli's  equa¬ 
tion  at  the  ocean  surface 
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The  linearized  forms  of  the  dynamic  equations  may  be  ob¬ 
tained  by  separating  the  effects  of  the  infinitesimal  surface 
waves  and  surface  current.  We  do  this  by  writing  the  velocity 
potential  and  surface  displacement  as, 

*  *  Y  +  ♦  (2.5) 

and 

C  *  X  +  h  (2.6) 

where  Y  and  x  are  associated  with  the  current  and  $  and  h  are 
associated  with  the  infinitesimal  surface  waves.  In  the  approx¬ 
imation  we  are  considering,  the  effect  of  the  surface  waves  on 
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the  current  is  negligible.  Therefore,  H'  and  X  satisfy  the  free 
surface  equations  independent  of  the  surface  height  (h)  and 
velocity  potential  (<J>),  i.e., 

4»  +  i  yy . yy  +gX  =  0  at  z  =  X  (2.7) 

t  2 

and 

X.  +  V'P'VX  =  P  at  z  =  X  •  (2.8) 

t  S  S  0  z 


Note  that  Equations  (2.7)  and  (2.8)  refer  to  the  unperturbed 
ocean  surface,  that  is,  unperturbed  by  the  surface  waves. 


To  determine  the  interaction  between  the  gravity  waves 
and  current,  we  substitute  Equations  (2.3)  and  (2.6)  into 
Equation  (2.2)  to  obtain, 

$  +  i-  V<J)  •  V<|>  +  V^-VY  +  gh  =  |  V4- » VT  +  gx)  at  z  -  X+h 

(2.9) 


The  right  hand  side  of  Equation  (2.9)  would  apparently  vanish 
using  Equation  (2.7)  except  that  the  two  expressions  are  not 
evaluated  at  the  same  boundary.  It  is  necessary,  therefore,  to 
expand  Equation  (2.9)  about  the  ambient  surface,  i.e.,  the  sur¬ 
rounding  region  on  which  there  are  no  gravity  waves.  The  ex¬ 
pansions  for  the  velocity  potential  are  as  follows: 


VY 


+  h 


z  =  X+h 


9 

9  z 


(VY) 


z  =  X 


(2.10) 
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ana 


« 

X+h 


X 


(2.1 


in  which  we  have  kept  only  terms  linear  in  the  surface  wave 
height  (h) . 


We  define  the  current  as  follows: 


VY  =  U  +  V  k  (2.12) 

s 

A 

where  k  is  a  unit  vector  in  the  vertical  direction.  The 
horizontal  current  U  is  significant  regardless  of  the 
generating  mechanism,  the  vertical  component  V  ,  however, 

5 

is  not.  For  example,  if  the  current  is  produced  by  a  non- 
dispersive  internal  wave,  then  |u|  >>  jv  |  so  that  only 
the  horizontal  component  need  be  retained  in  Equation  (2.12). 
If,  however,  the  current  is  produced  by  a  passing  swell,  the 
vertical  component  of  the  current  may  be  comparable  to  the 
horizontal,  so  that  both  must  be  considered.  The  particular 
generating  mechanism  will  therefore  determine  the  form  of  the 
current  used  in  the  analysis. 

Using  the  definition  of  the  current  (Equation  (2.12)]  and 


the  expansions  given  by  Equations  (2.10)  and  (2.11)  in  the  in¬ 
teraction  equation  (Equation  (2.9)]  we  obtain, 


<|>t  +  -j  V 4> •  V<t>  + 


g  + 


3V 

c 

W 


+  ( V cj)  +  V4»)  *vv 


h  +  W  •  V<*>  =  0 


at  z  =  x 
(2.13) 


If  we  now  restrict  our  analysis  to  the  linear  case,  i.e., 
linear  in  both  (j>  and  h,  we  reduce  Equation  (2.13)  to 

/  dV  \ 

<j>  +  Vf  *V(j>  +  ( g  +  )  h  =  0  at  z  =  X 


(2.14 


We  note  that  in  addition  to  the  standard  linear  interaction 
term  between  the  surface  waves  and  current  ( *  V 4> )  there  is  a 
second  interaction  due  to  the  vertical  component  of  the  current 
field.  This  modification  of  the  gravitational  acceleration 
has  been  discussed  recently  by  M.  Milder"^. 


In  determining  the  linear  form  of  the  kinematic  boundary 
condition  [Equation  (2.3)],  we  encounter  the  same  difficulty 
as  with  the  interaction  equation,  i.e.,  going  from  the  perturbed 
to  the  ambient  surface.  We  expand  the  terms  in  Equation  (2.3), 
just  as  before,  and  obtain 


V 


s 


z  =  X+h 


9z  ) 


(2.1:0 


to  terms  linear  in  h.  Using  Equation  (2.15)  along  with  the  other 
expansions  [Equations  (2.10)  and  (2.11)]  and  the  kinematic  bound¬ 
ary  condition  for  the  current  [Equation  (2.8)]  in  Equation  (2.3) 
yields , 
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(ft +  7# 


,7  +  74».y  +  7V  •  VJ h 
s  s  s  s  s  gf 


v$*vx  -  vv  *vx 

s  S  S  s  S 


at  z  =  X 


which  to  terms  linear  in  <J>  and  h  becomes 


dh 

art 


3V 


+  h  -  V*'7X  at  z  =  X 


(2.16) 


s  s 


The  incompressibility  condition  applied  to  the  surface  current 
yields  the  relation 


9  V  3  2m,  2  -*• 

S  _  g  V  =  _n  ‘w  =  -7  .  U 

3z  "  az2 


which  when  substituted  into  Equation  (2.16)  results  in, 

(ft  +,s5)h  ‘  If  -  7#'7SX  at  *'x  •  (2’171 


We  may  replace  the  gradient  of  the  ambient  surface  in  Equa¬ 
tion  (2.17)  by  taking  the  gradient  of  Equation  (2.7)  and  using 
the  definition  of  the  current  to  write 


VX 

s 


■  -  5  »sk + 15(32  + 

*=  -  |  f  +  V  7  V  1 

g  I  dt  s  s  sj 

Substituting  Equation  (2.18)  into  Equation  (2.17)  yields, 
(ft +7s-“)h  “  +  ? 


(2.18) 


(2.19) 
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Neglecting  terms  quadratic  and  higher,  except  when  they  involve 
a  coupling  of  the  surface  waves  to  the  surface  current,  we 
obtain  from  Equations  (2.14)  and  (2  19), 


h 


t 


t 


+  VY‘V<t>  +  gh  =  0 

+  V  t  •  V  h  +  (V  •  U )  h  = 
s  s  s 


(2.20a) 


(2.20b) 


evaluated  at  z  =  X.  Although  these  equations  have  the  same 
structure  in  one  dimension  as  those  derived  by  Zachariasen  (1972) 
they  are  evaluated  at  z=x  and  not  z=0  as  were  his  expressions. 

The  small  terms  obtained  by  further  expanding  the  set  (2.20) 
about  z=0  are  found  to  be  of  importance. 

Milder  (1972) 6  has  shown  that  the  simplicity  of  these 
equations  can  be  maintained  by  noting  that  the  right-hand  side 
of  Equation  (2.20b)  can  be  written  as, 

il  ix.  ii  =  ii  ,, 

3z  3x  3x  3x  3 n  ) 


in  a  single  horizontal  dimension  x,  where  s  is  the  distance 

along  the  surface  z  =  y  (x,t)  and  is  the  derivative  normal 

A  dn 

to  that  surface.  Milder  then  shows  that  because  the  velocity 
potential  satisfies  the  Laplacian,  in  two  dimensions, 


V2« 


/  3 .  3 \  /  3_ 

'  3x  1  3z  /  V  3x 


0 


(2.22) 


that 


9  s  9  <j>  _  _  .  9  ft 
9x  9n  *  1  9x 


(2.23) 


and  Equation  (2.20)  reduces  to, 

<ft  +  gh  =  0 

(2.24) 

it  h  +  uxh  -  -1  +x 

where  the  x  derivatives  are  constrained  to  the  surface 

d  d  9  3 

z  =  x(x,t)  and  ~  is  the  Lagrangian  derivative,  =  9t  +  U  9x * 

This  analysis  may  be  extended  to  two  dimensions,  resulting  in 

9  TJ 

the  loss  of  the  term  in  Equation  (2.20b)  (see  Appendix  in 
Watson  and  West,  PD-73-048) . 
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3.  MODAL  ANALYSIS  WITH  A  HARMONIC  SURFACE  CURRENT 

This  section  is  concerned  with  the  interaction  between 
surface  gravity  waves  and  a  current  distribution 

U  =  U  (x  ,  t)  ex 


where 

U (x, t)  =  S  UK  cos  K (x-Cjt)  (3-1) 

K 

and  e  is  a  unit  vector  in  the  direction  of  propagation.  The 
approach  used  is  that  developed  in  Section  I  in  which  eigen- 
node  equations  were  constructed  to  represent  the  interaction 
of  gravity  wa/es  with  gravity  waves.  We  will  extend  that 
analysis  to  include  the  interaction  of  surface  gravity  waves 
to  surface  currents  prescribed  by  Equation  (3.1).  The  linear 
solution  to  the  gravity  wave-sur  Cace  current  interaction 
problem  has  been  obtained  by  a  number  of  investigators  using 
a  variety  of  approaches;  e.g.,  Rosenbluth  (1971),  Ko  (1971) 
and  Zachariasen  (1972).  In  this  section  we  will  indicate  in 
the  context  of  the  eigenmode  model  how  these  various  calcula¬ 
tions  differ  in  their  level  of  approximation.  In  the  next 
section,  we  will  extend  the  calculation  into  the  nonlinear 
regime . 

The  first  order  equations  for  the  velocity  potential 
and  surface  elevation  derived  in  the  preceding  sections 
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[Equations  (2.20)  and  (2.21))  easily  accommodates  the  current 
distribution  defined  by  Equation  (3.1).  At  the  unperturbed 
surface,  the  interaction  may  be  written  in  two  dimensions  as 


♦t  +  gh 


0  cos  k£ 
X 


and 


cos  K5  - 


Kh  sinK£ 


(3.2) 


where  5  *  x  -  Cjt  and  Cj  is  the  phase  velocity  of  the  internal 

wave  packet.  Because  the  velocity  potential  satisfies  Laplace's 
equation,  we  may  write  the  two-dimensional  Fourier  transforms 
of  <{>  and  the  surface  elevation  as 


<P  (r ,  z ,  t) 


<J>£(t)  d2k 


and  (3.3) 

h(r,t)  =  -t-  f  eilc'r  hf(t)  d2k 

^  L  k 

and  r  =  (x,y).  Although  the  wavenumbers  are  continuous  in 
Equation  (3.3),  we  will  find  that  the  dynamic  equations  restrict 
their  values  to  a  discrete  set.  The  resulting  expressions  will, 
therefore,  conform  to  the  definitions  introduced  in  Section  I. 

The  interaction  equations  may  be  expressed  in  terms  of 
the  mode  amplitudes  as 
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Lr 


+  gh->  +  i|k|  U  $-  cos  k£  e  d  k  =  0 


K 


K  k 


and 


Jao 


I 


h->  -  |  k  |  e  '  k  1  z  4u 

k  k 


+  i|*l  S 


K 


L'  h  +  cos  K£ 
K  k 


Eil-r  2 

K  U  h*  sin  KE,  e  d  2 


ik •  r  ,2, 
e  d  k 


(3.4) 


We  recall  that  the  linearized  equations  have  been  obtained  at 
the  surface  z  =  X-  The  exponential  terms  in  Equation  (3.4) 
ray  therefore  be  see  equal  to  unity  since  | k [ x  c<  1  for  the 
surface  gravity  waves  of  interest.  If  we  multiply  Equation  (3.4) 
by  exp  [-  *k'rj  and  integrate  over  the  surface  plane  (r) ,  we 
obtain 


K  +  9hj 


-Si 


K 


2  UK 


-iKcTt  iKcTt 

(kx  -  K)  e  +  (kx  +  KM’  J+K  e 


and 


h-  -  |k|  <p£ 


(3.5a) 


El  -iKcTt 
r  •  1 


<VKI  H-i  +Khj-S 


+  e 


iKc_t 

i 


(kx  +  K)  hWi  ~  KhM 


as  the  linearized  mode  coupled  equations  to  second  order  in  the 


interaction. 


3A.  The  Eigenmode  Expansion 

We  now  introduce  the  eigenmode  variables  constructed 
in  Section  I: 

b(i)(*>  -  A  (k)  |*j;  *  i  ~ 1  (3.6) 

where  u>(k)  is  given  by  the  linear  dispersion  relation  (- 
and  A(k)  is  chosen  to  be 


A  (k) 


|g|2 

/2tT  o>(k) 


(3.7) 


By  taking  the  sum  and  difference  of  Equations  (3.5),  we  obtain 
the  eigenmode  equations 


ib{±)(£)  ?  uj (k) b  {±)  (jc)  =  A (k)  jiF.(jt)  *  **1*1  r  (jt) \ 

\  1*1  2  ) 

(3.8) 

where  F^k)  and  F2(k)  are  defined  by  the  right-hand  sides  of 
Equations  (3.5a)  and  (3,5b),  respectively. 

We  may  rewrite  Equation  (3.8)  in  the  more  explicit  form 


3b  (*}  (it) 
5t 


+  ioi(k)b(+)  (it)  * 


...K.  .  +  «>(k) 

w (k-K) 


K _ m(k) 

kx~K  W (k-K) 


k-K  -iKc_t 
e  1 

|  it-K  | 


(£-k) 


(S-K) 
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1  - 


k  +K 
x 


u  (k) 

u (k+K) 


ib(-)  (k+K) 


(3.9a) 


and  the  second  coupled  equation 


ib(“} (k) 
dt 


-  i  u(k)  b1*1  (it)  =  -  i  |it 


T-  K 


k-K  -iKcTt 


-►  ~r 

k-K 


1  +  ^  >  b (-)  (ic-K) 


kx-K  +  ,0  (k-K) 


1  + 


k  -K 
x 


<*>  (k) 
a)  (k-K) 


(  +  )  -*• 
hy  ’ (k-K) 


i  y  tI  vk  iKcifc 

4  °K  |k|  6 

K  I  k+K  | 


1  - 


K 


4-  +  (k) 

k.  +K  oj  (k+K) 


b(_) (k-K) 


1  - 


k  +K 
x 


u)  (k) 
w(  k+K) 


bK  1  (k+K) 


(3.9b) 
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where  we  have  used  Q  =  c^K  as  the  frequency  of  the  internal 

wave.  Examining  the  structure  of  the  coupling  coefficients 

in  Equation  (3.9)  makes  it  clear  that  and  b^  ^  are 

only  weakly  coupled,  since  |£|,  kx  >>  K  .  It 

is  only  necessary,  therefore,  to  discuss  the  equations  for 

either  b^  or  b^  but  not  both. 

We  may  remove  the  explicit  time  dependence  from 

* 

Equation  (3.9)  by  introducing  the  variables 


c(s)  <K>  -  b<‘>  (S)  eiVlt 


(3.10) 


If  we  introduce  Equation  (3.10)  into  Equation  (3.9a)  and 
use  the  decoupling  approximation,  we  obtain  the  equation 
for  right  travelling  waves 


ic(1c)  *  [a)  (k)  -  k^j}  c(£) 


VK 


S^-.i  +  *  +  J-jSL  c(S-i 

c-K|  kx  K  w(k'K)  , 

L 


_K _  w(k) 

k+K  co  (k+K) 


(k+K) 


where  we  have  deleted  the  extraneous  (  +  )  superscript.  A 
second  and  perhaps  stronger  justification  for  the  form  of 
Equation  (3.11)  is  based  on  the  notion  of  a  resonant  inter¬ 
action  between  the  surface  gravity  waves.  This  notion  was 

This  transformation  was  suggested  by  Kenneth  Case  and  the 
resulting  Hermitian  equations  are  explored  elsewhere,  Watson 
West  and  Case,  1973  (III)  . 


o  delate  urms  whose  frequencies  could  not  be 


thiri  or«.Ier  interaction.  A  similar  argu 


ca*.  vanisn  con 


Equations  *2.20)  and  (2.21),  the  dynamic  equation  may  be  added 
directly  to  the  nonlinear  expression  developed  in  I. 

The  composite  equation  in  one  dimension  is 


(k-K)  +  a.  ,„c(k+K)] 


upling  coefficients  arc  given  in  i.  the 

..n  ir  Equation  (3.11)  and  we  have  neglected  wind 


3B.  Approximate  Two-Dimensional  Equation 

First  let  us  note  that,  while  ic  is  a  continuous 
variable,  Equation  (3.11)  is  really  a  discrete  set.  Thus 
for  a  given  c(Jc)  ,  we  need  only  consider  the  equations  for 
c(ic')  where  it*  )c  +  ^  m^  K  with  m^  being  a  positive  or 
negative  integer.  To  obtain  a  Hermitian  form  of  Equation 
(3.11),  we  restrict  our  attention  to  the  case  of  a  single 
internal  wave  of  wavenumber  K;  then 


f  (ic) 

UK 

lie  1 

kx'K 

■  r 

lKl 

|*-*l 

g  (it) 

UK 

lie  1 

kx+K 

li  -  — 

"  T 

lkl 

k-K  | 

r  vk 

<d  00  \ 

srrarj 


and  introduce 
p  (ic)  - 


ui(k)  -  kx  Cj 


(3.13) 


(3.14) 


We  consider  the  condition  where  |£|  >>  K  which 

is  applicable  to  the  ocean  situation.  In  this  limit,  Equation 

(3.11)  reduces  to 


i  3c (k) 

at 


[w  (k)  -  k  c  ]  c  (ic)  + 

X  X 


| c (k+K) 


+  c(k-K) 


(3.15) 


The  change  of  variable 
c  (k)  =  t|,(ic) 


in  Equation  (3.11>  could  also  be  used  to  obtain 


i  3lKk) 

0  t 


=  [u  (k )  -  kxcI]  ij>(k) 


♦  £  £  L 

K  * 


(k+K)  j  (k+K) 


+  a^  (k)  y(k-K) 
ft  i 


(3.16) 


where 


aK(k)  =  j/|k|  |k  -  K  |  .  ( 3 . 17 ) 

which  reduces  to  Equation  (3.15)  for  |k|  >>  |k,. 

Equation  (3.16)  gives  the  interaction  of  the  gravity  waves 
with  an  internal  wave  spectrum  in  a  Hermitian  representation. 

To  obtain  an  alternate  expression  for  the  two- 
dimensional  problem,  we  focus  attention  on  Equation  (3.15) 
for  the  case  of  a  single  internal  wave.  We  note  that  we 
have  been  concerned  with  the  wavenumber  region  kj  K, 
so  that  we  may  employ  the  following  approximation 
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|  c(&+K)  +  c(k-K)j  nt  c  (it)  +  j  K2  Vj  c(£)  +  ...  (3.18) 


where  V+  is  the  gradient  operator  in  £-space.  If  in  addition 
to  Equation  (3.18)  we  look  at  the  eigenvalue  problem 


c(ic)  =  4»,(ic)  e"lXt 


Equation  (3.15)  becomes 


2  -►  2 
v;  (k)  + 

k  X  K2U, 


h»*  (k)  -  1 

1*1 


4>x  (k) 


where 


w'  (k)  =  a)  (k)  +  |)c  |  U.  -  k  c 


x  I 


(3.19) 


(3.20) 


(3.21) 


The  quantity  w'(k)  is  the  frequency  of  the  surface  wave  in 
the  coordinate  system  translating  with  the  phase  velocity  of 
the  internal  wave  and  is  therefore  a  constant. 

The  efficacy  of  using  the  expansion  given  by  Equation 
(3.18)  is  discussed  elsewhere  (III)  where  an  exact  one-dimensional 
eigenvalue  problem  is  solved  and  compared  with  that  obtained 
using  the  expansion.  The  eigenvalues  and  eigenfunctions  of 
the  approximate  expression  are  found  to  be  quite  close  to 
those  given  by  the  exact  solutions.  We  are  therefore  confident 
that  the  general  behavior  of  the  solution  to  Equation  (3.15) 
may  be  determined  by  examining  the  solutions  to  Equation  (3.20). 
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3C.  Approximate  One-Dimensional  Equation 

Let  us  from  here  on  focus  our  attention  on  the  one¬ 
dimensional  form  of  Equation  (3.15).  In  one  dimension  k  -  (>V0 
and  since  the  interaction  equations  form  a  discrete  set 
we  use  the  wavenumber  kx  =  nK,  n  =  0,1,  ...N-l,  d ,  N  +  l, 

For  a  single  internal  wave  of  wavenumber  K,  we  may  rewrite 

Equation  (3.15)  as 


1C 


n 


loj(n)  -  n/Jc 


n 


?VK 


'n+l 


+  c 


n-l| 


(3.22] 


where  f)(=KcI)  is  the  frequency  of  the  internal  wave.  We 
recall  that  the  strongly  interacting  surface  modes  are  those 
for  which  the  frequency  resonance  condition  given  by  Equation 
(3.12)  holds.  If  we  assume  that  the  central  surface  mode  is 
located  by  an  integer  n  >>  1  and  w(n)  n,  then  we  may  find 
a  simple  approximate  expression  for  Equation  (3.22). 

(i)  Expansion  in  Terms  of  Block  Waves 

We  introduce  a  new  variable  which  is  shifted  in  phase 

from  the  central  mode  n  =  N 

i[w-NR]t  (3.23) 

a  =  c  e 


into  Equation  (3.22)  to  obtain 

[gj  -  w(n)  -  Nft  +  nft]av  =  j  U0K  n  ^av  +  l  +  av-l^ 


22 


(3.24) 


Note  that  the  new  expansion  coefficients  [a^]  are  constants 
in  time  and  v  measures  the  deviation  of  the  wavenumber  from 
the  central  mode,  i.e.,  n  =  N  +  v. 

We  now  proceed  to  expand  u(N+v)  about  tu(N)  and  deter¬ 
mine  the  value  of  ft  in  terms  of  N.  We  select  the  frequency 
cj  to  be 

oj  =  uj(N)  +  6o)  (3.25) 


so  that  Equation  (3.24)  becomes 

L  +  i  (^)2  /5HkJ  av  -  ^  (N+V)  (av+1  +  av.l) 


or  expanding 


(N+v)  A, 


we  have 


av+l  +  av-l 


(3.26) 


where  ft  has  been  set  equal  to  ZgK/N  .  We  note  that 
Equation  (3.26)  is  the  expression  obtained  by  Rosenbluth  (1971) 
for  the  interaction  of  discrete  modes  when  the  expansion  set 
for  the  surface  elevation  and  velocity  potential  are  Block 
Waves.  It  might  also  be  mentioned  that  the  expansion  given 
by  Equation  (3.18)  in  the  two-dimensional  case  is  equivalent 
to  adding  zero  to  Equation  (3.26)  in  the  form  2av  -  2av  and 
taking  the  limit  of  continuous  v  as  done  by  Rosenbluth.  We 
also  note  that  although  terms  to  order  (v/N)  have  been  main¬ 
tained  in  Equation  (3.26),  terms  of  lower  order  have  previously 
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been  deleted,  i.e.,  order  v/N.  These  are  the  terms 


r- : — —  in  the  coupling  coefficients  [Equation  (3.13)] 

x 

which  arise  from  the  current  gradient  terrs  in  Equation  (2.21) 
and  were  first  discussed  by  Zachariasen  (1972). 

Because  of  the  inconsistency  in  the  order  of  terms 
maintained  in  Equation  (3.26),  we  take  a  slightly  different 
approach  than  that  of  Rosenbluth.  We  introduce  the  variable 


av(t) 


cn(t) 


i  [w  (N)  -Nil]  t 


1V2 


(3.27) 


into  the  one-dimensional  form  of  Equation  (3.11)  to  obtain 
av(t)  =  -  i  [  cju  ( n )  -  w(N)  -  v.<]  av(t) 

\  n-l  oi  (n-1 )  J  v-1 

-  (l - iy  +  fy?-)  a  .  (3.28) 

\  n+1  w  (n+1) J  v  +  1 

If  we  expand  oj(n+v)  about  w(N)  and  equate  with  j  v  gK/N 
then 

jj(N+v)  53  0,  (N)  +  vft  ;  n  =  N  +  v  (3.29a) 


and 


w  ( n )  _ .  i  4.  1 

wTntTT  ~  2n 


(3.29b) 
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Substituting  Equation  (3.29)  into  Equation  (3.28)  results  in 
the  expression 


av(t) 


U  K 
o 


n  (a 


v-1 


-  av+l’ 


3U  K 

+  — —  —  (a  ,  +  a  , ) 
8  v  v-1  v+1 ' 


(3.30) 


which  is  exact  to  order  .  The  second  term  in  Equation  (3.30) 

'  1 

is  of  order  smaller  than  the  first  and  invites  a  perturbation 
approach  to  the  problem. 

(ii)  Perturbation  Solution  of  Equation  (3,30) 

Let  us  consider  a  solution  to  Equation  (3.30)  of  the 

form 

av(t)  =  a^o) (t)  +  6ay(t)  (3.31) 

where  6a  (t)  is  of  order  ^  smaller  than  a/°^  (t) .  if  we  scale 
v  N  v 

the  time  in  Equation  (3.30)  as 

T  »  U  KNt  (3.32) 

o 

and  introduce  the  perturbation  expression  Equation  (3.31), 
we  obtain 

37av0)(t)  -  7  [av-i(t)  '  av°l(t)]  (3'33a) 
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and 


where  we  have  deleted  terms  of  order  1/N2.  Equations  (3.33a 
and  b)  are  correct  to  order  1/N. 

It  is  clear  that  Equation  (3.33a)  is  the  recursion 
relation  for  Bessel  functions.  The  solution  to  our  first 
order  equation  can  therefore  be  written  as 

av°)(t)  *  £  (3.34) 

£=-□0 

where  Jm(i)  is  the  Bessel  function  of  the  first  kind  of  order 
m.  The  expansion  coefficients  (d^o))are  determined  by  the 
fact  that  only  m  =  0  Bessel  functions  are  non-zero  at  time 
x=0  and  are  selected  to  satisfy  the  initial  conditions  for 
the  problem. 

Using  the  first  order  solution  given  by  Equation  (3.34) 
and  the  additional  recursion  relation  for  Bessel  functions 


2n 

T 


J  (t) 
n 


Jn-1<T> 


Jn+l(t) 


in  Equation  (3.33b)  we  obtain 


(3.35) 
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00 

a  , 
dr  Sav 

1  1 
J  i 

6av-l  "  6av+l 

£=-00 

*Vv  l-v 

V  —3 -  +  . .  J 

dx  T 


£-v 
(3.36) 


Equation  (3.36)  is  clearly  a  linear  equation  in  6ay  with  a 
driving  force  given  by  the  sum  over  Bessel  functions.  We 
may  solve  Equation  (3.36)  by  rewriting  the  equation  as 


3T  «VT>  '  J  [6av-l(T)  '  6Vl(T)]  *  Fv(T) 


(3.37) 


where 


VT)  ’  4N  ",  ~l 


L  a<°> 


aj,.v(t)  ,.v 

v  -ar — +  V 


(3.38) 


If  we  introduce  the  Laplace  transforms 

<£[6av(x)]  =  6av(s)  =  f  e“ST  6av(x)dx 


and 


(3.39) 


^[Fv(x)J  =  Fvis)  =  J  e_ST  Fv(x)dT 


into  Equation  (3.37),  we  obtain 

1 


s  6av(s)  -  f  I6av-1(s)  -  6av+1(s)]  *  Fy(s)  +  6ay(o) 


after  integrating^ I -T-rMby  parts 


1d6a  I 


(3.40) 
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We  observe  that  Equation  (3.40)  may  be  integrated 
by  using  the  discrete  Green's  function  g(v,s)  which  satisfies 
the  equation 


sg(v,s)  -  t  [g(v-l,s)  -  g(v+l,  s)] 


vo 


(3.41) 


The  solution  to  Equation  (3.41)  is  given  by 


g  (v,  s) 


-  £  f 


-ivKx  , 
e _ dx 

s-i  sin  Kx 


(3.42) 


which  may  be  verified  by  direct  substitut  ton .  The  formal 
solution  to  Equation  (3.40)  may  now  be  written  as 

6a  (s)  =  S  g(v-v',  s)  I" F v /  ( s)  +  6av#(o)l  .  (3.43) 

V  \T  /  ^ 


The  inverse  of  6av(s)  may  be  obtained  by  recalling 

X"1  j^[f1(T)]X[f2(T)]J  =  J  f  x  (t ' )  f  2  (x-'C ' )  dx' 


(3.44) 


where  -1{f}  is  the  inverse  Laplace  transform  of  f.  Using 
Equation  (3.44)  on  Equation  (3.43)  yields 

■vi:  g  (v-v#  ,  x-x,)[Fv  (x')  +  6av/(o)  6  (x  1 )  ]  dx  ' 


6av(x) 


(3.45) 
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The  Green's  function  in  Equation  (3.45)  is  given  by 


g (v-v '  ,  i-  ') 


L  j  -  i(v-v  ' )  Kx  ^i  :  ')  sinKx 
/Trf  -«> 


(3.46) 


Because  we  are  dealing  with  a  discrete  set  of  states,  we 
approximate  the  limits  on  this  integral  by  ±  which  is  the 
periodic  condition  imposed  on  the  surface  waves  over  the 
wavelength  of  the  internal  wave.  In  this  approximation, 
Equation  (3.46)  yields 

g  (v-v' ,t-t ')  «s  Jy_v,iT-  ')  (3.47) 


so  that  Equation  ( 3 . 43 )  becomes 


00 


l-v' 


X,- V  ' 


Vv'  (T~c')dT' 


as  the  general  solution  to  Equation  (3.33b)  with  given 

by  the  initial  conditions  for  the  solution  to  (3.33a). 

The  surface  elevation  is  given  by  the  discrete  form 
of  the  one-dimensional  Fourier  transform  in  Equation  (3.3). 
This,  together  with  the  defining  equations  ((3.6)  and  (3.10)] 
yields 
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h  (x ,  t) 


ik  (x-cTt) 
n  I 


ik  (x-cTt) 
n  I  ' 

e 


(3.49) 

We  note  that  for  a  situation  with  only  right  travelling  waves, 

i.e.,  predominantly  in  the  direction  of  the  wind,  ^  0 

for  k  >0  and  c  ^  ^  =  (c^)*so  that 
n  -n  n 


h  (x ,  t) 


-  z  > 


vr-  sin  [k  (x-c  t)  +  6  ] 
„  pi  ^  II  x  n 

n->0  n 


(3.50) 


where 


c(+) (t) 
n 


i0 


r  e 
n 


n 


(3.51) 


Using  Equations  (3.27),  (3.34)  and  (3.48),  we  have  for  the 

slope  variable 


c  (t ) 
n 


(t)  +  6ay(t) 


(io(N) 


NI2)  t 
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(3.52) 


where  we  have  set  <5a  (o)  =  0.  Inserting  Equation  (3.52) 
into  the  expansion  for  the  surface  elevation  results  in 
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h  (x ,  t) 
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(3.53) 


where 


(3.54) 


=  vK (x  -  Cjt) 


Equation  (3.53)  is  a  superposition  of  travelling  wave  with 
time  dependent  amplitudes. 

If  we  expand  the  sine  terms  in  Equation  (3.53),  we 
find  the  surface  elevation  can  be  written  as 


h (x , t)  =  | G (x , t) |  cos  [0N  -  arc  tan  (Gj/Gr) ]  (3.55) 


where 


G  (x ,  t )  = 


E  -fc  E 

N+V>  0  £=-oo 


<N+*)  + 


»  je  i 


g(v-v  ,t-t  )  Fv,(t')  di'  x  exp  i(kv£  -  v  j) 


(3.56) 


The  function  G(x,t),  therefore,  determines  the  modulation 
of  the  central  mode  k^(  =  NK)  .  We  may  express  the  envelope 

function  in  terms  of  cn(t)  rather  than  the  surface  elevation 
so  that 


G(x,t) 


CN+v(t)  itV  +  N  t) 

- -  e 


N+v>0  N+v 


(3.57 


If  we  use  the  slope  variables  introduced  in  T,  i.e., 


(  +  )  c(  +  )  i  [w  (n)  -n.2]  t 
n  n 


(3.58) 


then  we  would  have  from  Equation  (3.57) 
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where 


(3.59) 


cg  =  (oj  (N+v)  -  ui  (N)  ]/k 


(3.60) 


is  the  group  velocity  of  the  surface  wave  packet. 

The  modulation  function  for  the  surface  slope 
may  be  found  in  a  similar  manner  to  be 
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If  we  introduce  the  variable 


i(KC+j) 

z  *  e 

into  Equation  (3.61)/  then  we  have  an  expansion  in  terms 
of  zv.  This  result  can  be  used  for  formal  manipulation  of 
the  interaction  equation  when  the  solution  is  expressed  as 
a  generating  function, 

av  *V 
v 
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4.  NUMERICAL  CALCULATIONS 


4A.  Introduction 

In  this  section  we  wish  to  explore  the  relative 
effects  on  the  redistribution  of  energy  within  a  given 
initial  spectrum  by  both  the  mechanism  of  nonlinear  inter¬ 
action  of  surface  gravity  waves  and  the  interaction  of 
these  waves  with  a  prescribed  surface  current.  To  do  this 
we  use  the  full  nonlinear  expression  shown  in  Section  3, 

Equation  (3.13).  Our  discussion  will  cenv  r  on  two  calculations, 
both  giving  some  insight  into  the  relative  importance  of 
the  interactive  mechanisms.  The  first  calculation  will  be 
that  of  an  impulsively  applied  harmonic  surface  current 
interacting  with  a  single  surface  mode.  The  growth  of  the 
spectrum  as  a  function  of  time  will  be  discussed,  as  well  as 
the  nonlinear  effects  on  the  surface  modulation.  This  sim¬ 
plified  picture  will  be  compared  with  the  linear  calculations 
of  the  Lewis  experiment  as  an  indication  of  the  veracity  of 
the  numerical  techniques  used  in  the  integration.  The  second 
calculation  replaces  the  single  surface  mode  of  the  preceding 
case  with  a  discrete  surface  spectrum  interacting  with  an 
internal  wave  with  typical  oceanographic  parameters. 


4b.  Linear  Surface  Modes 


In  this  calculation,  we  wish  to  focus  on  the  inter¬ 
action  between  a  prescribed  surface  current  and  surface 
qravity  waves.  To  isolate  this  interaction  from  the  nonlirear 
interaction  between  surface  gravity  waves,  we  consider  small 
amplitude  surface  waves,  so  that  the  nonlinear  interaction 
is  inhibited.  The  analysis  in  Section  3  is  appropriate  for 
calculations  in  which  the  parameters  of  the  problem  are 
characteristic  of  the  ocean  environment,  it  is  inappropriate, 
however,  for  conditions  found  in  typical  tank  experiments. 

This  situation  becomes  clear  if  we  examine  the  relative  size 
of  the  interaction  coefficients  to  the  doppler  shift  in  the 
i genmode  frequency. 

We  recall  that  the  frequency  of  the  (N+v)  sideband 
was  approximated  in  Section  3  by 

">(N+v)  «  uj(n)  +  vfi  .  (4>1) 

I'be  next  order  term  in  this  expansion,  which  was  neglected 
in  Section  3,  would  be  v2ft/4N,  where  ^  =  w(N)/2N  and  the 
central  wavenumber  is  given  by  kN  =  NK.  Typical  values  for 
ocean  and  tank  experiments  are 
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Ocean 


Tank 


n  -  10-2  sec-1 

K  -  2*10-4  cm-1 

N  ~  500 

Uq  -  2  cm/sec 

v2  </4N  -  5V  10~6  v2 

U  k  ~  0.1  sec-1 
o 


^  1  sec 
K  ~  0.05  cm  1 
N  ~  8 

-2 

U  ~  10  cm/sec 
o 

V2  -V4N  -  0.03  v2 

U  k  —  2><10  ^  see 
o 


In  the  ocean  environment  the  interaction  coefficient  is 

4 

typically  10  larger  than  the  ioppler  shift.?  howoy«r, 
in  tanks  the  doppier  shift  is  typically  in  order  of 
magnitude  larger  than  the  coupling  coefficient.  The  tank 
is  therefore  in  a  separate  region  of  the  experimental  range 
than  is  the  ocean. 

To  see  this  result  in  detail,  we  consider  the  inter¬ 
action  equation  for  the  surface  gravity  waves  with  a  surface 
current  generated  by  a  single  internal  wave  with  the  frequency 
correction  term  included 


a  (t) 
v 


=  IV 


4N 


av(t) 


1  U  N  K 

2  o 


v-1 


(t)  - 


av+l (t) 


(4.2) 

The  terms  neglected  in  Equation  (4.2)  are  an  order  of  magnitude 
smaller  than  the  coupxing  coefficient  which  in  turn  are  an 
order  of  magnitude  smaller  than  the  frequency  correction  term. 
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Let  us  consider  the  problem  of  an  impulsively-applied 
harmonic  current  to  a  single  small  amplitude  surface  wave, 
i.e.,  at  time  t  =  0  we  "turn  on"  the  interaction  between  a 
single  surface  mode  and  internal  mode  both  of  infinite  spatial 
extent.  If  we  concern  ourselves  with  only  the  growth  of  the 
first  two  side  bands,  we  obtain  the  coupled  system 
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(4.3) 


where  k  =  NK  is  the  wavenumber  of  the  initial  wave.  We 
N 

may  combine  the  equations  in  (4.3)  'to  obtain 
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(4.4) 
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-5 

The  coefficient  of  the  forcing  function  is  of  the  order  10  sec 
and  the  solution  to  the  system  may  be  obtained  to  a  good 
approximation  b,  sf Lting  the  right-hand  side  equax  to  zero, 
especially  for  the  linear  case,  yielding 

a*l'fc)  =  “=J—  e  8N  sinot;  a  =  §  ^(^KUq)2  +  (jjj) 


when  a  t ,  lo)  =  0 . 

The  modulation  function  for  the  surface  elevation 
is  given  by  Equation  (J.57)  which  in  terms  of  the  av ’ s 
becomes 


(4.5) 


G(x,t)  = 
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/=-! 


(-i) 


-  v  ivK£  pi<  N 
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(4.6) 


where 


=  V  "  “(N)t 


and 


(4.7) 


£  =  x  -  crt 

x 


In  terms  of  the  solutions  for  the  impulsively  applied 
internal  wave ,  we  have 


a  (o)  el(’N  I  NKU 
u;x,t)  =  -2— -  Jl  - 
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1  “iTrTt 

°  siin  at  I  sin  K£  +  cos  K£  e 


[* 


(4.8) 


fo**  the  modulation  of  the  initial  wave.  For  the  early  time 
doluiii  i,  i.e.,  t  1,  we  can  write 

-  cos  Y  \ 


G  (x  ,  t) 


h  (o )  e 
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KU  t 
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sin  K4  + 


since 

write 
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(4.9) 


0.02 


-] 

sec 


so 


that 


t  <<  50  sec. 


We  can  also 


|  G  (x  ,  t ) 


%  h (o) o*  N 
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U  Kt 
o 


j^sin  K •' 


Or  CO£: 


(4.10) 

which  agrees  with  Zachariasen ' s  calculation  tor  the  initial 
value  problem.  Equation  (4.10)  may  now  be  used  to  check  the 
numerical  calculation  at  early  times. 

For  the  numerical  calculation,  we  selected  the  wave- 
number  of  the  surface  mode  to  be  =  0.412  cm-1  and  internal 
mode  number  K  =  0.0515  cm  ^ .  These  values  wore  chosen  to  be 
in  the  range  used  by  Lewis  (1971)  in  his  interaction  experi¬ 
ment.  For  an  internal  wave  of  phase  speed  of  Cj  -  24.39  cm/sec, 
the  surface  and  internal  wave  are  in  resonance.  Figure  (1) 
depicts  tne  discrete  spectrum  of  surface  waves  at  five  different 
times;  the  "0  sec"  graph  illustrating  the  initial  conditions 
for  our  problem.  The  initial  slope  |q  |  =  0.0524  is  that 
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of  the  mechanically  generated  surface  wave  in  the  Lewis 
experiment.  The  spacing  between  modes  is  given  by  the 
internal  mode  number  Ak  -  K. 

The  spectral  "snapshots"  in  Figure  (1)  indicate  the 
persistence  of  a  rather  i arrow  spectrum.  Only  the  sidebands 
oi  the  primary  wave  have  gained  appreciably  in  energy  after 
a  time  interval  corresponding  to  twelve  cycles  of  the  central 
mode.  The  other  modes  seem  to  acquire  and  lose  energy,  but 
not  to  steadily  increase.  This  is  seen  more  clearly  in  the 
following  figure. 

In  Figure  (2)  we  indicate  the  continuous  growth  of  the 
spectral  mode  slopes  in  time.  It  is  clear  that  three  of  the 
nine  modes  in  the  calculation  dominate  the  characteristics 
of  the  spectrum.  In  hindsight  then,  the  inclusion  of  the 
other  six  modes  was  superfluous  and  their  removal  from  the 
calculation  should  not  materially  alter  any  conclusions  made. 
The  effect  of  the  internal  wave  on  the  surface  spectrum  is, 
therefore,  primarily  local  in  k-space,  its  width  being  approxi 
matcly  2K  where  K  is  the  internal  mode  wavenumber. 

In  linear  theory  for  the  case  of  resonance  between  the 
primary  surface  wave  and  internal  wave,  the  surf  ice  modulation 
should  be  stationary  with  respect  to  the  interna]  wave.  We 
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q  =  MODULUS  OF  SURFACE  SLOPE 


IOH 


Figure  2. 


may  characterize  the  interaction  as  Lewis  (1971)  did  in  his 
experiment,  by  defining  a  quantity  called  the  surface  dis¬ 
tortion  (S)  , 


S 


(4.11) 


which  using  the  linear  calculation  in  this  section  reduces 
to 


S  =  Uo  Kt  sinK:  +  -j"  cosKr ] 


(4.12) 


Because  we  have  a  resonance  condition,  a  given  segment  of 
the  surface  wave  interacts  continuously  with  a  given  phase 
point  of  the  internal  wave  as  they  both  propagate  down  the 
tank.  The  length  of  time  of  the  interaction  can  there¬ 
fore  be  measured  as  a  function  of  distance  down  the  tank, 
i.e.,  t  =  x/Cj.  Keeping  this  in  mind  and  recalling 
Zachariasen ' s  calculation  of  the  maximum  surface  distortion, 
we  rewrite  Eq.  (4.12)  as, 


(4.13) 


This  expression  differs  in  the  initial  slope  from  the  one 
given  by  Zachariasen,  but  as  he  pointed  out  the  problem  he 


considered  in  detail  was  the  boundary  value  problem  where¬ 
as  Eq.  (4.13)  corresponds  to  the  initial  value  problem. 

Eq.  (4.13)  does  agree  with  the  linear  calculation  of  Ko  (1971) 
and  also  the  initial  value  problem  of  SUchariasen,  the 
solution  of  which  he  left  in  general  form. 

In  Figure  (3)  we  compare  the  results  of  the  numerical 
calculation  of  the  interaction  with  the  linear  expression  in 
Eq.  (4.13)  and  the  results  of  the  Lewis  experiment  (1971). 

The  linear  theory  agrees  with  the  numerical  calculation 
for  a  surprisingly  lonq  time.  One  should  recall,  nowevf-r , 
that  the  parameters  only  restrict  the  interaction  time  to  be 
less  than  50  seconds  for  the  linear  case.  Near  this  fin.o* 
the  two  results  diverge,  the  linear  case  growing  quadradi- 
cally  and  the  nonlinear  numerical  calculation  saturating 
to  a  surface  distortion  of  approximately  0.65.  It  is  un¬ 
fortunate  that  for  the  experimental  parameters  selected 
that  the  tank  was  not  a  few  feet  longer  so  as  to  distinguish 
experimentally  between  these  two  calculations. 

In  the  light  of  these  results  one  could  reproduce 
in  the  present  context  the  complete  analysis  giver,  by  Ko 
(1971).  The  third  curve  in  Figure  (3)  is  an  example  of 
such  a  calculation  in  which  the  initial  wave  is  not  in 
resonance  with  the  internal  wave  and  therefore  the  surface 
distortion  is  inhibited.  It  is  not  the  similarities  with 
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LEWIS  EXPERIMENT 


DISTANCE  (feet) 


linear  theory,  however,  but  the  differences  which  are  of 

interest  here.  In  linear  theory,  the  resonant  wave  and 

those  close  to  resonance  are  unique  in  that  the  surface 

distortion  grows  without  limit,  whereas  the  nonresonant 

waves  lead  to  finite  surface  distortions.  In  general, 

however,  if  we  start  with  a  resonant  system  the  nonlinear 

interactions  will  detune  the  system  thereby  inducing 

saturation.  The  saturation  level  for  the  initially  resonant 

system,  will,  however,  be  maximum. 

In  Figure  (4)  the  effect  of  including  the  four  wave 

coupling  terms  in  Eq.  (3.13)  can  be  seen.  The  first  effect 

noted  is  the  increased  rate  of  growth  in  each  mode,  except 

of  course  the  primary  mode  which  has  an  increased  rate  of 

attentuation.  To  understand  this  effect  we  recall  the 

discussion  in  I  about  the  preferential  amplification  of 

the  Benjamin-Feir  sidebands.  The  sideband  frequencies 

which  had  maximum  amplification  were  uj.  =  (1+q^  ) 

oF  o  o 

where  q^  is  the  slope  of  the  central  wave  in  our  notation, 
o 

The  corresponding  wavenumbers  are  k^  =  kQ  (l  +  2q^_  )  which 

o 

for  the  initial  conditions  in  our  problem  are  k-^  =  k  (1+0.104), 

BF  o 

i.e.,  +.043  cm  1  from  k  .  It  is  coincidental  that  the 

o 

Ben jamin-Feir  sidebands  lie  very  close  to  the  sidebands  coupled 
by  the  internal  wave  to  the  primary  mode,  i.e.,  +.0515  cm  1 
from k  .  This  near  coincidence  of  the  wavenumber,  however, 
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accounts  for  the  marked  increase  in  the  qrowth  1  ate  of 


the  sidebands. 

The  second  effect  apparent  from  comparing  Figure  4) 
with  Figure  (3) is  the  loss  of  oscillation  in  the  mode  ampli¬ 
tudes.  This  occurs  because  of  the  chanqe  in  the  phase 
velocity  induced  by  the  nonlinear  terms,  i.e.,  a  nonlinear 
dispersion  relation.  To  see  this  effect  in  mode  detail 
and  to  determine  if  a  finer  grid  for  th°  wavenumbers  would 
significantly  chanqe  the  spectral  growth  indicated  in 
Figure  (4),  we  introduced  modes  at  k^^-K  and  *  jK  in  the 
present  problem. 

In  Figure (5)  the  growth  of  the  nine  modes  in  the 
latter  calculation  is  indicated.  The  initial  amplitudes 
for  this  calculation  are  slightly  different  from  the  pro¬ 
ceeding  but  this  difference  should  be  of  no  consequence. 

Of  importance  to  the  early  time  behavior  is  the  initial 
choice  of  phases  for  the  modes.  We  see  that  with  the  phases 
selected  the  sidebands  at  kQ+ 2K  initially  give  up  energy 
to  the  system.  This  trend  is  soon  reversed,  however,  so 
that  at  late  times  ( t  ^  2  5  sec)  the  growth  of  the  modes 
directly  coupled  to  the  internal  wave  is  not  noticeably 
affected  by  the  initial  choice  of  phases.  We  also  see  that 
the  modes  at  half-odd  multiples  of  K  which  only  grow  by 
means  of  the  nonlinear  interaction,  play  no  part  in  the 
spectral  development. 


Comparing  Figures  (3),  (4)  and  (5)  wt-  see  that: 

(i)  the  energy  transferred  directly  from  the  internal  wave 
to  the  surface  spectrum  is  localized  in  k-space,  (ii)  the 
four  wave  interaction  spreads  this  energy  out,  i.e., 
broadens  the  •.•ffect  of  (i),  and  (iii)  the  malor  effect 
of  the  four  wave  interactions  is  between  waves  directly 
coupled  to  the  primary  wive  via  the  internal  wave  due  to 
the  coincin-_Tc  '  of  the  Ben jamin-Feir  spacing  (.043  cm-*) 
and  the  wavor.ui.b jr  spacing  (.0515  cm-*). 
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4C.  Ocean  Spectral  Development 

In  this  calculation  we  wish  to  contrast  the  spectral 
development  in  an  r  c- •  .  environment  with  that  just  calculated 
for  a  tank.  We  consider  an  impulsively-applied  interaction 
between  an  i:.t  Cl  i  A  in  the  ocean  and  a  surface  spectrum. 

We  assume  the  internal  wave  to  be  '10^  meters  long;  i.e., 

-4  -1  -2  -1 

K  =  2><10  cm  ,  witn  frequency  ft  =  10  sec  so  that  its 

phase  velocity  is  =  50  cm/sec.  The  central  wavenumber  of 

the  surface  spectrum  is  assumed  to  be  in  resonance  with  the 

internal  wave  and  tnercfore  has  a  wavelength  of  64.18  cm; 

_•>  _i 

i.e.,  kQ ~  9.79  10  "  cm  .  The  amplitude  of  the  central  mode 
is  1  cm,  so  that  its  slope  (q^)  is  0.1. 

The  estimate  in  Section  4B  indicates  that  the  coupling 
of  the  central  node  to  the  rest  of  the  spectrum  in  the  ocean 
is  much  greater  than  in  the  tank  case.  This  implies  that 
the  energy  will  drain  much  faster  from  the  central  mode.  For 
this  reason  the  present  calculation  is  done  with  21  modes  so 
as  to  avoid  any  difficulty  with  the  finite  size  of  the  spectrum 
in  determining  the  characteristics  of  the  spectral  development. 
To  simulate  the  ocean  environment,  we  assume  the  remaining 
twenty  modes  to  have  an  initial  amplitude  of  10  i.e., 

two  orders  of  magnitude  smaller  than  the  primary  and  to  have 
random  phases. 
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Due  to  the  large  number  of  modes  in  this  calculation, 
the  representation  of  the  results  is  slightly  different  than 
in  the  tank  case.  If  we  omit  the  four  wave  coupling,  the 
spectrum  is  shown  in  Figures (6)  and  (7)  at  different  times. 
Although  the  number  of  modes  is  discrete,  in  Figure  (6)  we 
join  the  mode  amplitudes  at  a  given  time  to  aid  the  eye  in 
determining  the  spectrum.  The  "0  sec"  spectrum  is  the  initial 
state  of  the  problem  and  is  depicted  as  a  very  nai row  spectrum. 
After  the  interaction  is  "turned  on"  the  spectrum  is  seen  to 
broaden  and  develop  a  bi-modal  character. 

Figure  (7 )  depicts  a  continuation  ol  the  spectral 
cievelopment  initiated  in  Figure  (6).  It  can  be  seen  that  the 
energy  is  being  symmetrically  transferred  to  both  higher  and 
lower  wavenumbers  as  the  interaction  persists.  It  is  evident 
that  the  more  central  mode  amplitudes  are  oscillating  in  time 
and  as  the  bulk  of  the  energy  is  transferred  to  the  spectral 
extremes  the  energy  in  the  central  part  of  the  spectrum  is 
decaying.  Indeed,  if  we  look  at  the  linear  part  of  the  solu¬ 
tion  to  this  problem  as  given  by  Equation  (3.52),  we  have 

qn(t)  ~  qN(o)  JN-vtr) 

Then  using  the  asymptotic  form  of  the  Bessel  function,  the 
central  part  of  the  spectrum  decays  as  l//t  and  oscillates 
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with  a  period  T  =  r— jj~  which  is  approximately  5  seconds 

ouo 

in  our  problem.  This  is  the  structure  that  is  beginning 
to  appear  in  Figure  (7).  The  central  mode,  whose  frequency 
is  =  9.79  rad/sec,  has  a  period  of  0.641  sec  and  is 
depleted  of  energy  in  approximately  eighteen  cycles.  A 
somewhat  lesser  amount  of  energy  is  then  returned  to  the 
mode  in  a  somewhat  shorter  time. 

Figure  (8)  depicts  a  further  contrast  with  the  tank 
experiment.  In  this  figure  we  show  the  surface  distortion 
of  the  ocean  for  the  21-mode  calculation.  If  we  examine 
the  expression  for  the  modulation  function  of  the  surface 
slope  given  by  Eq .  (3.61),  retaining  only  the  linear  terms 

ao  have 


Gs (x,t) 


A 

53  qN(0)  VT)  Z" 

v=-A 


(4.14) 


i (Ka+j) 

where  z  =  e  ,  2A  is  the  spectral  width  and  N  >  A  >>  1. 

We  may  approximate  the  sum  in  Eq .  (4.14)  by, 


A 

2  Z  Jv  <T)  »  2  (z  '  t)j 

v=-A 


which  is  the  generating  function  for  Bessel  functions. 
Using  the  definition  of  z  and  Eq.  (4.15)  we  obtain  for  the 
slope  modulation  function, 
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Uo  =  2.0  cm. /sec. 
K  -2%  IO‘4/cm. 


TIME  (seconds) 


Gg  (x ,  t )  qN(o)  exp^ii  cosKi;]  .  (4.16) 

Because  Eq .  (4.16)  has  modulus  unity,  it  is  clear  that  to 

first  order  both  the  surface  distortion  and  slope  distortion 
are  constant.  This  linear  result  is  in  agreement  with  the 
numerical  calculation  indicated  in  Figure  (ft) 

For  the  initial  conditions  of  this  problem,  the  four 
wave  interaction  terms  do  not  affect  the  qualitative  results 
shown.  In  fact,  for  the  first  14  seconds  or  approximately 
twenty  periods  of  the  primary  wave  there  is  essentially  no 
quantitative  difference  in  the  surface  distortion.  Further 
calculations  using  more  "physical"  initial  conditions  will 
be  conducted  elsewhere. 
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4D.  Conclusions 


It  is  clear  from  the  calculations  discussed  in 
Section  4B  and  4C  that  the  ocean  and  tank  experiments  fall 
into  two  quite  different  interaction  regimes.  This  is  most 
clearly  illustrated  by  the  growth  of  the  surface  distortion 
found  in  the  tank  experiment  and  calculation  depicted  in 
Figure  (4)  and  the  complete  lack  of  growth  found  under 
oceanographic  conditions  for  an  equivalent  interaction  de¬ 
picted  in  Figure  (8 ) .  This  result  points  up  the  extreme 
caution  required  in  the  design  of  tank  experiments  whose 
purpose  is  to  induce  information  about  Interactions  on  the 
real  ocean. 

The  calculations  clearly  demonstrate  the  need  for 
detailed  analysis  of  experiments  to  determine  their  appli¬ 
cability  to  the  ocean  environment.  To  this  end  a  program 
of  calculations  is  being  initiated  which  will  determine 
spectral  modifications  produced  by  surface  currents  under 
a  variety  of  oceanographic  conditions.  Calculations  will 
also  be  made  for  the  corresponding  tank  experiments  where 
such  experiments  can  be  conducted.  An  example  of  such  a 
calculation  would  be  that  of  a  stationary  current  pattern 
in  the  lab  which  simulates  the  convergence  zone  of  a  long 
internal  wave.  Since  such  a  simulated  internal  wave  could 
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be  constructed  by  contouring  the  bottom  of  a  tank  appro¬ 
priately,  the  difficulty  encountered  in  generating  a  trans¬ 
lating  internal  wave  in  a  tank  is  avoided. 
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